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AEROSOL  EXTINCTION  TRANSMISSOMETRY  AT  VISIBLE  AND  INFRARED  WAVELENGTHS 


APPLIED  TO  SPECIALIZED  PROBLEMS  IN  CB  DEFENSE 


1.  INTRODUCTION  AND  BACKGROUND 

With  the  increasing  availability  of  high-speed  electronic  comput¬ 
ers  over  the  past  two  decades,  complex  calculations  of  optical  extinction 
and  scattering  of  electromagnetic  radiation  by  atmospheric  aerosols  using 
the  Mie  theory1  have  become  commonplace. 

The  technology  has  advanced  so  swiftly  that  specialists  measuring, 
e.g.,  cloud  extinction  of  aerosols  often  are  unaware  of  advances  being  made 
simultaneously  by  other  specialists  measuring,  e.g.,  angular  scattering 
patterns  of  aerosols  with  nephelometers.  Such  measurements  and  calcula¬ 
tions  can  be  done  in  numerous  wavelength  regions,  but  the  technology  which 
concentrated  first  upon  the  visible  wavelengths  now  is  expanding  rapidly 
especially  into  the  regions  of  relatively  high  atmospheric  transparency 
called  "windows"  in  the  infrared. 

This  technological  expansion  can  be  measured,  for  example,  by  a 
sampling  of  papers  appearing  one  journal  in  one  year  alone,  discussing  aer¬ 
osol  measurements  or  calculations  of  extinction, 2-12  angular  light  scatter¬ 
ing  or  nephelometry,8* 13-18  1 idar, 1 9—22  ancj  related  aerosol  technology 
including  thermal  emissivity23  and  calibration  of  optical  particle 
counters.24 

The  author  has  noted  several  trends  in  all  of  this  work.  First, 
the  utility  and  information  content  of  simple  extinction  measurements  has 
tended  to  be  overlooked  in  the  rush  to  exploit  scattering  techniques. 
Second,  the  polarization  techniques  do  not  seem  to  include  applicable 
techniques,  some  new,  that  can  be  applied  both  to  extinction  and  to  polari¬ 
zation  measurements.  Third,  that  if  taken  in  combination  at  this  point  in 
time,  the  overall  technology  might  be  further  advanced  than  is  realized  by 
workers  in  the  individual  disciplines. 

For  example,  remote  characterizat ion  or  quantitative  analysis  of 
tropospheric  aerosols  seems  a  distinct  possibility,  both  by  conventional 
transmission  (extinction)  spectroscopy  and  by  angular  scattering  or 
backscattering  techniques.  This  report  examines  these  possibilities  by 
first  reviewing  the  technology  of  aerosol  extinction  spectroscopy. 

The  Beer-lambert  equation  for  atmospheric  aerosols  can  be  written: 


-  In  TX  =  axCL, 


(1) 
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where  for  a  given  wavelength  (pm),  T\  is  the  fractional  transmittance, 
is  the  optical  mass  extinction  coefficient  of  the  aerosol  (m2g*1),  C  is  the 
aerosol  mass  concentration  (g  m-3),  and  L  is  the  optical  path  length  (m). 

Some  form  of  Equation  (1)  has  Deen  used  by  spectroscopists  virtually  since 

the  birth  of  spectroscopy  and  the  equation  has  become  so  well  known  that 

the  power  of  simple  extinction  measurements  of  aerosols  often  is  not  fully 

appreciated. 

For  example,  using  suitable  multi-wavelength  techniques, 25,26 
extinction  measurements  can  contain  sufficient  information  to  allow  direct 
determination  of  the  mean  size  of  a  spherical  droplet  aerosol  distribution, 27 
or  even  remotely  to  characterize  chemical  reactions  in  cloud  droplets 
at  isosbestic  wavelengths, 28  without  the  use  of  polarized  radiation  or  other 
more  sophisticated  techniques. 

It  is  not  surprising,  however,  that  polarization  techniques  have 
received  heavy  emphasis  since  their  exploitation  in  the  early  1940s  to 
characterize  military  fog-oil  smoke  mean  droplet  diameters  based,  in  part, 
in  earlier  work  at  the  U.S.  Bureau  of  Standards  by  Stratton  and  Houghton2^ 
using  the  Mie  theory^  that  was  first  put  forth  in  1908  but  had  to  await  the 
availability  of  high-speed  computers  for  exploitation. 

Stratton  and  Houghton  measured  the  polarization  ratio  only  for 
transparent  spheres  for  which,  if  the  complex  index  of  refraction  m^  is 
given  by  {n  -  ik)x  where  n*  is  the  real  index  and  kx  is  the  imaginary  index 
(closely  related  to  the  absorption  coefficient),  then  k^  =  0. 

Polarization  measurements  can  be  used  to  determine  particle  si2e 
distribution  in  addition  to  mean  size  of  aerosol  droplet  distributions. 

This  was  recognized  and  the  technique  was  developed  in  the  1960s.  By  1969, 
Kerker  and  his  co-workers30-32  had  extended  the  techniques  of  particle  size 
distribution  measurements  to  the  extent  discussed  in  his  classical  book  on 
light  scattering33  that  gives  an  excellent  review  of  the  technology  as  it 
existed  chen. 

In  recent  years,  there  is  an  increased  awareness  that  aerosol  par¬ 
ticle  absorption,  e.g.,  in  the  infrared  where  is  not  zero,  can  be  inter¬ 
preted  in  new  ways  using  not  only  extinction  measurements2^'2?  or 
polarization  measurements, 33  but  the  two  together. 26  the  author  concurs 
that  this  is  a  promising  field  for  new  research. 

2.  AEROSOL  EXTINCTION  MEASUREMENTS 

It  is  known  that  the  variation  of  from  values  near  zero  in  the 
visible  wavelengths  to  significant  ones  in  the  infrared  wavelengths  pro¬ 
vides  a  aethod  for  the  determination  of  aerosol  mean  droplet  size  a  mass 
concentration  from  simple  extinction  measurements  if  two  or  more  wave¬ 
lengths  are  selected  properly  for  observation. 

For  example,  at  the  He;Ne  laser  wavelength  X  =  0.63  pm,  kx  is  near 
zero  and  typical  water  fog  droplets  have  diameters  0^  such  that  DU>>A 
or,  to  use  Kerker ‘s  criterion,  the  size  parameter  xOy/X  >.  2.0  so  that  the 
conditions  for  geometric  scattering  exist. 


This  leads  to  a  constancy  of  the  product  cix-D^  — ►  3. 0-3. 2,  where 
a\  is  defined  in  Equation  (1).  At  the  same  time,  the  product  ax-Dp  decid¬ 
edly  is  not  constant  with  Dm  at  wavelengths  where  kx  is  not  equal  to  zero 
in  the  infrared.  This  is  illustrated  in  Figure  1,  calculated  from  the  Mie 
theory. 


In  Figure  1,  curves  are  shown  for  the  infrared  wavelengths  of  8.5, 

10.5  and  12.57  {or  4 it)  pm,  as  well  as  for  a  composite  of  wavelengths  in  the 
8.5-12.57  pm  region  where  extinction  can  be  measured,  e.g.,  by  a  simple 
broadband  transmissometer  operating  at  8-13  pm  in  the  atmospheric  window 
region  there. 

While  operation  at  a  specific  wavelength  such  as  the  10.6  pm  CO2 
laser  line  for  comparison  to  a  visible  wavelength  such  as  the  He:Ne  laser 
line  at  X  =  0.63  pm  has  advantages,  a  broadband  infrared  reference  band 
also  can  be  used.  In  fact,  related  findings  using  a  filtered  light  source 
in  the  visible  at  X  =  0.515  pm  and  a  9-12  pm  broadband  transmissometer  to 
study  the  formation  and  dissipation  of  water  fogs  were  reported  at  least  as 
early  as  1970  by  Carlon.34 

When  the  techniques  of  Carlon,  et  al,25,27  are  extended  to  use  the 
Mie  theory  to  calculate  functions  of  the  transmi ttances  T^  at  X  =  0.63  pm 
and  the  broadband  or  composite  wavelength  band  8.5-12.57  pm  vs.  water  drop¬ 
let  diameter  D^,  curves  like  those  shown  in  Figure  2  are  obtained.  It  is 
seen  that  the  curves  are  widely  separated. 

Thus,  a  sensitive  method  is  provided  for  the  determination  of  mean 
droplet  diameter.  Experimental  points  for  a  steam-generated  water  fog  are 
shown  in  Figure  2;  these  give  excellent  agreement  with  calculated  curves 
for  example  at  =  8  pm.  This  result  was  cited  oy  Weinman,  et  al,35  as 
approximating  Oeirmendjian’s  1969  Model  C-1  (or  vice-versa),3®  where  Op  is 
taken  as  the  mass  median  diameter. 

What  is  uncertain  is  the  kind  of  distribution  to  which  this  diame¬ 
ter  is  most  applicable.  Kerker33  has  given  an  informative  discussion  of 
particle  sue  distributions,  and  of  transmission  (extinction)  and  turbidity 
measurements.  In  this  report,  turbidity  can  be  taken  as  the  product  ox-C 
in  Equation  (1),  having  the  units  m'1. 

Justificati  for  the  use  of  a  composite  of  wavelengths  (Figure  1 
and  abscissa.  Figure  2)  can  be  found  by  examining  calculated  curves  and 
steam-generated  fog  data  points  for  the  individual  wavelengths  X  =  83, 

10.5  anti  12.57  (or  4*)  pm  and  plotting  curves  like  Figure  2  as  is  done  in 
Figure  3. 

Viewing  the  plots  from  top  to  bottom  in  Figure  3,  it  can  be  seen 
that  the  fog  data  gave  the  best  agreement  with  calculation  for  a  mass 
median  diameter  (MHO)  of  6  pm  at  X  3  8.5  p»,  HMD  *  10  u»  at  X  =  10.5  pm, 
and  NMD  =  6  p*  at  X  «-•  pa 

In  his  discussions  of  particle  sire  distribution  functions, 
kerker33  cautions  that  mean  particle  diameters  obtained  by  using  methods 
like  those  discussed  in  the  present  report  may  be  adequate  for  most  pur¬ 
poses  if  the  distribution  is  sufficiently  narrow,  but  strictly  are  valid 
only  for  monodisperse  systems  - 
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Figure  1.  Values  of  a^.Ojj  vs  0„  Calculated  Fro*  the  Mie  Progras  for 
Water;  the  “Constant"  Function  Actually  is 
Qaaped  Oscillation. 
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Figure  2.  Comparison  of  Computed  Curves  Relating  Visible  to  Infrared 
Optical  Transmittance  of  Water  Fog  to  Experimental  Data 
Points  for  Cooling,  Steam-Generated  Water  Fog  Clouds. 
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Figure  3.  See  Figure  2;  Curves  for  Three  Infrared  Wavelengths 
Comprising  the  Broadband  Composite  of  Wavelengths 
Shown  as  the  Abscissa  of  Figure  2. 


By  using  a  composite  of  wavelengths  to  obtain  a  mean  diameter,  one 
is  averaging  both  over  wavelength  intervals  and  over  particle  size 
distributions.  But  the  results  indicate  that  the  technique  works  very 
well,  especially  under  dynamic  droplet  growth  or  evaporation  conditions 
where  what  is  wanted  most  is  a  real-time  indication  of  approximate  droplet 
size,  e.g.  in  studies  of  developing  or  dissipating  water  fogs. 34, 37 

The  dependence  of  a\  upon  A  is,  in  fact,  highly  variable  oue  to 
variations  in  (n  -  ik)x,  in  addition  to  which  ct\  varies  with  Dj;  as  was  dis¬ 
cussed  previously.  This  is  illustrated  in  Mie  calculations  for  the  extinc¬ 
tion  coefficient  of  water  droplets  (a\,  ordinate)  in  Figure  4  vs  Dvi  for 
several  wavelengths.  At  certain  wavelengths,  e.g.,  A  =  12.5  pm,  tne  com¬ 
bined  contributions  of  refractive  (n\)  and  absorptive  (k\)  components  of 
water  droplet  extinction  lead  to  values  of  a\  that  virtually  are  independ¬ 
ent  of  droplet  size. 

Experiments^?  seem  to  confirm  for  example  that  ct-j 2.5  1S  nearly 
constant  with  Dp,  suggesting  that  the  liquid  water  content  C  in  Equation 
(1)  of  a  tropospheric  optical  path  could  be  monitored  by  a  simple  12.5  pm 
transmissometer  with  good  precision  tor  droplet  sizes  up  to  about  15  pm. 38 
In  this  application,  however,  care  must  be  taken  to  acccunt  for  absorption 
due  to  hydrogen-bonded  molecular  clusters  (water  clusters)  in  the  vapor 
leading  to  absorption  easily  confused  with  droplet  absorption. 5 

At  wavelengths  called  "isosbestic  points", 28  the  complex  indices 
(n  -  ik)x  of  some  liquids  are  such  that  a\  remains  constant  even  though  the 
liquid  solutions  comprising  the  droplet  aerosol  vary  widely  in  their  chemi¬ 
cal  composition,  e.g.,  for  droplets  of  orthophosphor ic  and  related  acids 
that  might  be  formed  by  burning  prosphorus  in  air  of  varying  relative 
humidity  (RH)  such  that  solute  (acid)  concentrations  vary  from  near  zero  to 
85%  by  weight. 

In  Figure  5,  A  =  11.4  pm  is  an  isosbestic  point  calculated  from 
refractive  index  measurements  of  pure  orthophosphoric  acid  (H3PO4)  in 
water. 39, 40  Real  atmospheric  smokes  obtained  by  burning  phosphorus  may 
contain  mixtures  of  acids  depending  on  rates  of  combustion  in  the  presence 
of  varying  humidities  and  so  may  not  give  experimental  spectra  like  those 
for  pure  H3PO4  in  Figure  5  where  an. 4  remains  essentially  constant  at 
0.13-0.15  m2g-i,  while  acid  concentration  varies  from  near  zero  to  85%  by 
weight.  This  will  be  investigated. 

Because  of  the  constancy  of  an. 4  compared  to  the  wide  ranges  of, 
say,  ag.7  (see  Figure  5)  with  solute  concentration,  remote  characterization 
of  the  aerosol  becomes  a  possibility  if  something  is  known  about  the  drop¬ 
let  constituents.  It  can  be  shown  that: 


ax  4nk\f  (n>x)/Ap, 


(2) 


where  p  (g  cm'3)  is  the  droplet  acid  (solution)  density,  i.e.,  mass 
density,  and: 
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Figure  4. 


Calculated  Curves  For  Water  Droplets  Showing  a\  vs  D^;  Note 
Constancy  of  ai2.5  at  All  Wavelengths  Up  to  About  =  15  pin. 
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Figure  5 


Computed  Values  of  a\  for  Liquid  Droplet  Aerosols  Compri 
Several  Concentrations  of  Orthophosphor ic  Acid  {H3PO4)  i 


sing 
n  Water. 
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(3) 


f(n»A)  = _ 

<nX2  +  kx2)2  +  4(nx2  -  kx2)  +  4 


where  mx  =  (n  -ik)x  and  the  function  of  Equation  3  is  plotted  in  Figure  6 
for  values  of  nx  marked  on  the  curves  and  ranging  from  0,1  to  4.0,  and  for 
values  of  kx  shown  on  the  abscissa.  It  is  emphasized  that  the  above  dis¬ 
cussion  applies  to  Rayleigh-scattering  smokes  where  D^«A  as  clearly  is 
the  case  for  orthophosphoric  or  other  acid  aerosols  (smokes)  produced  by 
burning  phosphorus  in  moist  air. 

Often  it  is  useful  to  plot  extinction  coefficients  or  fractional 
transmittances  vs  rather  than  vs  the  size  parameter  tiD^/A,  for  specific 
wavelengths.  Two  such  calculated  plots^l  are  shown  for  water  droplets  in 
Figures  7  and  8.  In  Figure  7,  the  curves  are  calculated  for  a  product 
C  x  L  =  5.0  (Equation  (1))  and  for  wavelengths,  reading  from  top  to  bottom 
of  the  left-hand  ends  of  the  curves,  of  A  =  0.63,  8.5,  10.5  and  4n  pm, 
i.e.,  the  same  wavelengths  considered  in  Figures  1-3.  In  Figure  8,  the 
curves  are  calculated  for  A  =  10  pm  to  show  the  functions  vs  of  total 
extinction,  and  of  the  scattered  and  absorbed  components  that  are  summed  to 
give  total  extinction,  i.e.. 


aT  =  as  +  aA. 


(4) 


where  the  subscripts  indicate  total  (T),  scattering  (S),  and  absorptive  (A) 
extinction. 

Sassen26  recognized  that  aerosol  cloud  extinction  measured  in  con¬ 
junction  with  angular  scattering  measurements  could  lead  to  remote  sensing 
of  cloud  composition  by  using  two  wavelengths  like  the  laser  wavelengths 
A  =  0.63  and  10.6  pm.  But  the  present  discussion  indicates  that  such 
remote  characterization  could  result  from  cloud  extinction  measurements 
alone  at  two  closely-spaced  infrared  wavelengths,  one  of  them  an  isosbestic 
point. 


Additional  data  might  be  obtained  by  combining  extinction  and 
angular  scattering  measurements  at  an  isosbestic  wavelength.  Optical 
parameters  can  carry  only  so  much  information,  but  the  idea  is  provocative 
enough  to  warrant  further  investigation. 

Other  investigations  of  remote  characterization  of  tropospheric 
aerosol  clouds  are  possible,  such  as  those  based  on  effects  peculiar  to 
Christiansen  wavelengths. 42  A  complete  discussion  is  given  by 
Carlon, 4*5,42  where  Reference  5  discusses  aerosols  ranging  in  size  from 
molecular  clusters  to  water  cloud  droplets. 


of  0. 1-4.0  Shown  on  the 
Abscissa,  ‘'on  Rayleigh- 
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Figure  7 


Fractional  Transmittance  for  a  C  x  L  Product  =50  (Equation  (1)) 
of  Water  Fog  vs  Oroplet  Oiameter  for  Wavelengths,  Reading  from 
Top  to  Bottom  of  the  Left-Hand  Ends  of  the  Curves,  of  A  =  0.63, 
8.5,  10.5  and  4n  pm.  Which  Are  the  Wavelengths  Also 
Considered  in  Figures  1-3. 
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Figure  8.  a\  vs  0^  Calculated  From  the  Hie  Program  at  A  =  10pm  for  Water 
Droplets;  the  Subscripts  T,  S  and  A  Refer  to  Ecuation  (4). 
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3. 


POLARIZATION  AND  ANGULAR  SCATTERING  MEASUREMENTS 


Kerker33  has  given  an  excellent  review  of  techniques  used  to 
determine  particle  size  distributions  using  measurements  like  these,  and 
the  present  author  will  not  elaborate  except  to  raise  a  few  pertinent 
points. 


Sassen26  verified  that  angular  scattering  patterns  calculated  from 
the  Mie  theory  can  be  observed  rather  closely  by  experiment  not  only  for 
spherical  particles,  but  for  irregularly-shaped  ones  as  well.  Thus  polari¬ 
zation  or  angular  scattering  techniques  combined  with  some  of  the  recent 
advances  in  extinction  measurement  technology  discussed  here  should  provide 
new  combinations  for  the  remote  characterization  of  tropospheric  aerosols. 

The  polarization  method33  for  evaluating  the  distribution  of 
sphere  sizes,  utilizes  polarization  of  scattered  radiation  obtained  from  a 
monochromatic  light  source  at  various  angles  of  observation.  The  intensity 
of  the  components  of  the  scattered  light  whose  electric  vector  vibrates 
perpendicular  and  parallel  to  the  plane  of  observation,  1 1 ( 0 )  and  1 2 ( 9 ) #  is 
measured  at  a  number  of  angles,  where  0  is  the  scattering  angle  measured 
from  the  direction  of  the  incident  light  beam.  The  polarization  ratio, 
P(0),  is  defined  such  that: 


P(8)  =  I2(9)/I 1 (9) • 


(5) 


For  example,  P(8)  can  be  calculated  as  a  function  of  0  for  any 
assumed  particle  distribution  such  as  the  zeroth  order  logarithmic  distri¬ 
bution  (ZOLD)  and  for  any  standard  deviation,  o0,  -anging  from  aQ  =  0 
(i.e.,  a  monodisperse  aerosol),  to  values  of  0.3  or  more  where  the  curves 
flatten  and  the  method  approaches  its  limit  of  usefulness. 

Figure  9  shows  the  family  of  curves  calculated  for  the  He:Cd  laser 
wavelength  A  =  0.4416  pm,  0^  =  0.3  pm,  n^  =  1.434  and  kx  =  0  so  that 
rax  =  (1,484  -  i.Q)x  and  the  density  is  p  =  0.98  g  cm-3.  The  curves  are  for 
dioctyl  phthalate  (OOP),  a  liquid  that  is  commonly  used  for  test  aerosols.6 
Agreement  between  calculation  and  measurement  using  the  polarization  method 
is  good  to  excellent,  often  within  a  few  percent. 

The  effect  of  the  wavelength  of  the  monochromatic  light  upon  the 
functions  shown  in  Figure  9  is  of  special  interest.  This  has  been 
investigated30*-*1  for  visible  and  near-infrared  wavelengths,  including 
cases  such  as  vanadium  pentoxide  (Vp0<>)  for  which  strong  absorption  occurs 
at  shorter  visible  wavelengths,  i.e.,  where  kx  is  not  zero.3^ 

From  the  discussion  in  the  present  report  of  special  techniques 
using  extinction  measurements  in  the  infrared  where  kx  is  not  zero,  it  is 
suggested  that  combinations  of  visible  and  infrared  wavelengths  can  opti- 
mi/e  desired  parameters  using  polarization  techniques.  Since  the  size 
parameter  nO./A  depends  on  wavelength,  a  change  in  the  wavelength  of  the 
incident  ligrit  will  cause  a  shift  in  the  position  of  the  curve  peaks  in 
figure  9  such  that  larger  values  of  A  shift  the  peaks  to  the  right. 
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Figure  9 


Polarization  Ratio  P (G)  vs  Scattering  Angle  0  for  OOP. 


2 


Multi -wavelength  observations  have  been  used  to  show  that  particle 
size  distributions  obtained  at  different  wavelengths  are  consistent, 33  but 
when  wavelength  spacing  is  as  great  as  X  =  0.63  and  10.5  pm  additional  work 
is  needed  to  confirm  this  and  also  to  determine  whether  the  useful  range  of 
a0  can  be  extended  beyond  0.3  in  measurements  using  the  polarization  method. 

The  scattering  ratio  method  for  the  determination  of  particle  size 
distribution  is  a  variation  where  the  angle  of  observation  is  held  constant 
and  the  polarization  ratio  is  measured  or  calculated  over  a  spectrum  or 
over  a  number  of  discete  wavelengths. 

The  method  was  developed^,  44  in  the  early  1960s  and  was  perfected 
by  using  monodisperse  polystyrene  latex  spheres  and  known  combinations  of 
single  sphere  sizes  to  give  mixtures  of  known  size  distributions  checked  by 
electron  microscopy. 

Thus  the  scattering  ratio  method  gives  an  alternative  way  of 
studying,  e.g.,  the  variation  of  a0  at  some  Dp  with  X  at  some  constant  9, 
although  the  same  thing  can  be  done  by  using  the  polarization  method 
directly  and  plotting,  e.g..  Figure  9  with  X  as  the  abscissa  at  constant  0. 

Turbidimetric  particle  size  distribution  methods  invariably  are 
similar  to  Wallach,  et  al.46  Nelson47  proposed  a  method  by  which  calcu¬ 
lated  scattering  cross  sections  from  the  Hie  theory  are  compared  to  experi¬ 
mental  extinction  data  using  monodispersed  polystyrene  latex  spheres  in 
distilled  water  as  a  calibration  means,  and  a  computer  program  is  used  to 
null  the  error  vs  wavelength  thus  forcing  a  unique  solution  for  the  size 
distribution  which  is  teken  as  bimodal  to  account  for  coagulation.  The 
method  also  is  said  to  yield  particle  volume  concentration  and  the  scatter¬ 
ing  cross  section  function,  Q^.  Extrapolation  of  the  technique  to 
tropospheric  aerosols  remains  to  be  demonstrated. 

4.  SUMMARY  AND  CONCLUSIONS 

This  paper  has  concentrated  on  advances  that  have  been  made  in  the 
remote  characterization  of  tropospheric  aerosols  by  extinction  and 
turbidimetric  methods  over  the  nearly  two  decades  that  have  elapsed  since 
the  publication  of  Kerker's  classic  text.33  The  technology  has  been 
updated  by  including  results  not  previously  published  with  improved  tech¬ 
niques  in  polarization  and  angular  scattering  methods.  Perhaps  most 
importantly,  the  author  has  tried  to  support  the  argument  that  if  all  meth¬ 
ods  presently  available  in  this  technology  are  taken  into  account  at  this 
time,  our  ability  to  remotely  characterize  or  quantitatively  analyze  atmos¬ 
pheric  aerosols  is  further  advanced  than  workers  in  the  individual  light 
measurement  disciplines  realize. 

It  is  concluded  that  by  combinations  of  well-established  and 
recently-developed  light  scattering  techniques,  investigators  can  undertake 
promising  new  areas  of  research.  Some  examples  include  combinations  of 
extinction  and  polarization  measurements  with  computer  programs  to  study 
isosbestic  or  Christiansen  wavelengths,  and  optimization  of  desired 
parameters  as  functions  of  wavelength  or  particle  size  distribution. 
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APPENDIX  A 


LASER  MONITORING  OF  MASS  CONCENTRATIONS 
OF  MONODISPERSE  TEST  AEROSOLS 


A1 .  INTRODUCTION 

Aerosol  sampling  usually  is  considered  to  be  a  very  complicated 
problemA1  in  aerosol  physics.  It  may  be  desired,  for  example,  to  determine 
the  efficiency,  as  a  function  of  particle  diameter  (pm),  with  which  aer¬ 
osol  particles  of  various  kinds  are  drawn  into  an  aspirated  sampling  tube 
or  head  or  into  the  human  respiratory  system. 

The  principal  experimental  difficulty  is  one  of  establishing  and 
maintaining  a  test  aerosol  of  known  characteristics,  including  composition, 
mass  or  number  concentration,  and  particle  size  distribution.  The  aerosol 
sampler  is  placed  in  the  test  aerosol  and  is  operated  at  a  known  flow  rate 
for  a  known  period  of  time,  after  which  a  determination  is  made,  usually  by 
weighing  a  filter  paper  in  the  sampler,  of  the  mass  of  the  collected 
aerosol.  This  is  ratioed  to  the  particle  mass  contained  in  the  sampled 
volume  of  the  test  aerosol  to  determine  the  sampler  efficiency. 

The  problem  of  maintaining  the  test  aerosol  is  not  a  trivial  one. 

It  must  be  sampled  periodically  to  determine  if  (and  how)  it  is  changing 
with  time  as,  for  example,  particles  settle  or  agglomerate.  This  sampling 
disturbs  and  draws  material  from  the  test  aerosol,  and  the  subsequent  anal¬ 
yses  usually  require  at  least  several  minutes  to  complete,  during  which 
period  the  characteristics  of  the  test  aerosol  may  continue  to  change. 

This  appendix  describes  a  technique  by  which  the  test  aerosol  mass 
concentration  is  continuously  sampled  by  a  laser  beam.  This  requires  that 
the  optical  constants,  particle  size  distribution,  and  other  properties  of 
the  aerosol  particles  be  known.  The  test  aerosol  materials  discussed  in 
this  paper  are  liquids  that  are  generated  as  spherical  droplets,  for  which 
accurate  (•*■  or  -  10%)  optical  calculations  can  be  made  using  the  Mie  theory. A2 

For  approximation  purposes  in  optical  calculations,  it  is  often 
useful  to  treat  an  aerosol  as  if  it  consisted  entirely  of  monodisperse  par¬ 
ticles  of  some  equivalent  diameter. A3  But  a  better  approach  is  to  generate 
monodisperse  droplet  aerosols  directly,  as  was  done  in  the  work  reported 
here. 
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A2 


THEORY 


In  the  geometric  optical  scattering  regime,  particle  diameters  are 
much  larger  than  the  wavelength  A  (pm)  of  observation.  That  is,  D^^A. 
Under  these  conditions,  it  is  straightforward  to  showA4  that: 


<*A  =  3Qx/2DMp 


(Al) 


where  a  (m2g-1)  is  the  mass  extinction  coefficient  of  the  spherical 
particles,  p  (g  cm-3)  is  the  particle  density,  and  Qx  is  the  crcss-section 
efficiency  factor,  which  is  very  nearly  constant  and  equal  to  2.0  for  the 
liquid  aerosol  droplets  discussed  in  this  paper, A3  so  that: 


ax 3/D Mp 


(A2) 


i he  mass  extinction  coefficient  is  used  in  the  Beer-Lambert 

equation: 


-  In  TA  =  axCL, 


(A3) 


where  Tx  is  the  optical  transmittance  at  laser  wavelength  A,  C  (g  m-3)  is 
the  aerosol  mass  concentration,  and  l  is  the  optical  path  length  uniformly 
filled  by  the  aerosol  (m).  C,  the  quantity  sought  for  the  test  aerosol,  is 
then: 

C  =  *•  — i_  In  Tx'*'-  S pfi  In  Tx  (A4) 

aAL  3  L 


Several  liquids  used  to  generate  test  aerosols  in  work  described 
in  this  appendix  are  shown  in  the  table. 


Tsbl*.  liqjokk  U**d  to  Ctonsrsto  TmI  AMOK* 


Symbol 

Substance 

P. 

density, 

(g/enP) 

Twnjwr- 

ature 

CO 

Keal 

iiuiox* 

<"o) 

DMH 

o-Dimethyl 

phihaiate 

U92 

20.8 

i.si&a 

DliP 

p-lhvthjl 

phthiis'.e 

s.sa 

17.7 

1.5029 

IMiP 

nDitnrtyl 

phihalate 

i.o*a 

25.0 

1.4925 

i)OI* 

Ihoctyl 

phthsiate 

<h&3 

20.0 

1.4550 

*  Tb«  sodium  D  line  t*  at  (ho  wavelength,  *r.<J  the  real 

wd«  «t  thf  Hp-Ne  User  waveienglh  nl0.63.um «*  nearly  identical  to 
this  value  fur  usftsfssreni  colorless  iujukl*. 
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Mie  calculations  were  performed  for  the  He:Ne  laser  (X  =  0.63  pm), 
with  results  as  shown  in  Figure  A1 ,  and  it  was  confirmed  that  the  approxi¬ 
mation  of  Equation  (A2)  is  quite  precise,  i.e.,  that  the  functional  depend¬ 
ence  of  ao.63/  the  extinction  coefficient,  on  Dy  is  very  -■''milar  for  all 
liquid  aerosol  materials  considered,  with  only  liquid  densities  p  having  a 
predictable  effect. 

The  laser  wavelength  (0.63  pm)  lies  very  close  to  the  sodium  D 
line  (0.589  pm),  and  the  materials  listed  in  Table  A1  are  transparent  and 
colorless  in  the  visible  wavelengths  (i.e.,  they  have  a  real  refractive 
index  n  but  a  negligible  absorption  coefficient  or  imaginary  index  k). 

Thus  it  was  very  convenient  to  obtain  the  real  indices  of  these  materials 
using  a  standard  laboratory  ref ractometer  at  the  sodium  D  line  and  to  use 
these  in  Mie  calculations  for  the  0.63-pm  wavelength.  The  error  introduced 
by  doing  this  is  very  small,  as  is  indicated  by  Figure  A2,  which  compares 
the  ax  vs  Dy  functions  for  DMP  calculated  at  X  =  0.589  and  0.63  pm. 

After  preliminary  testing,  it  was  found  convenient  to  use  the  liq¬ 
uid  DOP  (Figure  A3)  in  most  subsequent  work.  Since  p  =  0.98  g  cm~3  for 
DOP,  the  approximations  in  Equations  (A2)  and  (A4)  can  be  still  further 
simplified: 


a0.63'~'3/Dp* 


(A5) 


and: 


C/^-  9m-  1n  t0.63 
3  L 


(A6) 


It  can  be  seen  from  Figure  A3  that  the  approximation  of  Equation  (A5) 
is  quite  precise  in  the  geometric  regime  (Dy»X). 

A3.  EXPERIMENTAL  PROCEDURE 

A  1-m3  test  chamber  was  constructed  and  used  as  shown  schemati¬ 
cally  in  Figure  4.  The  He-.Ne  laser  was  optically  aligned  through  pinholes 
in  opposite  chamber  walls.  The  pinholes  were  utilized  to  prevent  window 
effects.  A  pressure-equalization  scheme  between  the  chamber  interior  and 
room  air  prevented  aerosol  particles  from  flowing  in  either  direction 
through  the  pinholes  during  testing. 

Although  Figure  4  shows  an  aerosol  cloud  coming  from  the 
generator,  the  chamber  was  stirred  during  testing  to  ensure  uniform  mixing. 
This  was  important  because  laser  path  length  L  of  the  laser  beam  in  the 
test  chamber  (which  was  1  m  long)  must  be  filled  uniformly  with  the  test 
aerosol  during  measurements. 

The  aerosol  generator  was  of  the  spinning  disk  design^  to  produce 
essentially  aonodisperse  liquid  droplets  of  test  aerosol.  The  droplet  size 
0y  disseminated  by  the  spinning  disk  generator  is  deterained  by  the  diame¬ 
ter  and  angular  velocity  of  the  disk,  by  the  liquid  flow  rate,  and  to  a 
lesser  extent  by  other  paraaeters  including  teaperature. 
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EXTINCT 


Figure  A1 .  Values  of  an  „  vs  Dn  Calculated  Fret  the  Mi©  Theory  tor 
Liquids  listfta  in  Table  Al. 
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Figure  A2.  Values  of  o\  vs  Ojj  Calculated  From  the  Hie  Theory  for  Dimethyl 

Phthalate  (DHP),  dsing  the  Sodium  D  Line  Refractive  Index  Measured 
at  Its  Wavelength  of  0.589  pm  (Dashed  Curve)  and  Assuming  the  Same 
Index  at  the  He: We  Laser  Wavelength  A  a  0.63  pm  (Solid  Curve). 
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Figure  A3.  Values  of  0053  vs  Calculated  Fro»  Hie  Theory  for 
Oioctyl  Phthalate  (OOP). 
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AEROSOL  PARTICLE  SIZING  BY  USER 

(AFTER  UNIFORM  MIXING) 


0  to  100 

PARTICLE  DIAMETER  (ytn) 

OR  (DISC  VELOCITY)-1 


Figure  A4.  Scheaatic  Representation  of  the  Experimental  Setup;  the  Chamber 
Was  Stirred  to  Ensure  Uniform  Aerosol  Mixing  in  the  Laser  Beaa. 
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Aerosol  droplet  size  was  checked  by  microscopic  examination  of 
samples  deposited  on  glass  slides.  The  generator  was  easily  capable  of 
disseminating  mass  concentrations  of  test  aerosols  in  the  range  C  =  0. 1-1.0 
g  m-3  for  these  tests.  Other  monodisperse  aerosol  generators  are  commer¬ 
cially  available,  although  some*  do  not  have  the  aerosol  output  capability 
of  the  unit  used  in  these  tests.  Plastic  microspheres**  also  were  used  in 
some  trials,  for  which  the  ag.63  vs  Op  functions  were  not  much  different  from 
those  shown  in  Figure  A1. 

For  testing,  aerosol  samplers  were  mounted  completely  inside  the 
chamber  with  their  inlets  positioned  perpendicular  to  and  directly  adjacent 
the  laser  beam.  This  ensured  that  the  laser  beam  and  inlets  intercepted 
identical  samples.  A  refinement  to  this  procedure  involves  using  a  second 
laser  beam  perpendicular  to  the  first  and  nearly  intersecting  it,  thus  pro¬ 
viding  laser  beam  "cross  hairs",  which  can  be  positioned  directly  in  front 
of  the  sampler  inlet  tube  under  test. 

A4.  RESULTS  AND  DISCUSSION 

The  purpose  of  this  appendix  is  to  describe  a  technique  by  which 
the  mass  concentration  of  a  sized  test  aerosol  can  be  continuously  sampled 
by  a  laser  beam  to  improve  experimental  procedures.  Only  a  brief  discus¬ 
sion  of  typical  test  data  for  aspirated  sampler  efficiencies  will  be  given 
here.  Figure  A5  presents  some  data  for  OOP  droplets  and  illustrates  the 
methodology  discussed  in  this  paper. 

In  Figure  A5,  the  solid  and  dashed  curves  labeled  Aspirated  Sampler 
ae  typical  of  the  efficiencies  measured.  For  these  aspirated  samplers,  the 
sampling  efficiency  falls  to  less  than  10%  in  the  DOP  droplet  size  range 
D(|  =  30-40  pm.  The  He:Ne  laser  can  be  considered  100%  efficient  in 
sampling  the  test  aerosol,  since  laser  attenuation  depends  entirely  on 
extinction  by  the  aerosol  droplets  (e^  =  100%). 

Multiple  scattering  of  the  laser  beam  is  not  a  problem,  since  the 
test  aerosol  concentration  is  always  chosen  to  give  high  transmittances 
(0. 7-1.0),  thus  conserving  the  aerosol  material.  The  assumption  that  the 
actual  mass  concentration  of  the  test  aerosol,  Cactual.  is  equal  to  that 
determined  from  the  laser  using  Equation  (A6),  Ciaser,  is  valid  to  within 
the  error  of  the  Hie  calculations,  which  was  found  to  be  as  good  as  plus  or 
minus  10%  in  this  application. 

The  equivalent  mass  concentration  collected  by  the  aspirated 
sampler,  CS3jBp'|eq,  is  determined  by  weighing  the  filter  paper  in  the  sam¬ 
pler  before  and  after  an  aerosol  trial  and  is  reproducible  within  about 
plus  or  minus  10%.  The  percent  efficiency  of  the  sampler  e$  is  then 
100  x  (CsaWpiecj/Cactuai).  When  the  test  aerosol  particle  diameter  is  only 


*  The  Berglund-Liu  monodisperse  aerosol  generator  is  marketed  by  Thermo- 
Systems,  Inc.  (TSI).  P.0.  Box  3394,  St.  Paul,  Minnesota  55165. 

**  Microspheres  of  glass,  polystyrene,  and  other  materials  are  available 
from  Ouke  Scientific  Corp.,  445  Sherman  Avenue.  Palo  Alto,  Calif.  94306. 
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DETERMINATION  OF  ASPIRATED  SAMPLER  EFFICIENCY 

USING  LASER  BEAM 


=  PARTICLE  DIAMETER,  pm 


Figure  AS.  Typical  Sampling  Efficiencies  for  COP  Test  Aerosols  Are  Shown 
by  the  Curves;  the  Methodology  of  Measurements  Also 
Is  Illustrated. 
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a  few  micrometers  or  less,  the  aspirated  sampler  efficiency  approaches 
100%,  and  these  samplers  can  be  used  directly  to  confirm  laser  measurements 
of  aerosol  mass  concentration  as  was  reported  previously  for  smoke 
aerosols.A6 

AS.  CONCLUSIONS 

The  discussion  in  this  appendix  has  been  limited  to  the  determination 
of  aspirated  aerosol  sampler  efficiencies  using  laser  reference  techniques. 
These  techniques,  of  course,  are  not  limited  to  sampler  efficiency  measure¬ 
ments  but  can  be  used  wherever  it  is  desirable  to  continuously  monitor  the 
mass  concentration  of  a  test  aerosol  of  spherical  droplets,  e.g.,  in  stud¬ 
ies  of  the  human  respiratory  system. 

By  using  a  He:Ne  laser  in  the  visible  wavelengths  (A  =  0.63  pm), 
aerosol  measurements  are  made  in  the  geometric  scattering  regime  (D^  A) 
using  materials  for  which  many  simplifications  and  precise  approximations 
are  possible  from  the  Mie  theory.  The  method  is  fast,  accurate,  and  elimi¬ 
nates  the  need  for  precalibration  of  standard  samplers  against  which,  for 
example,  candidate  samplers  must  be  compared. 
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APPENDIX  B 


NASS  EXTINCTION  COEFFICIENTS  ESTIMATED  FOR 
NONABSORBING  SPHERICAL  AEROSOL  PARTICLES 
IN  THE  GEOMETRIC  SCATTERING  REGIME 


The  Beer-Lambert  equation  for  aerosols  can  be  written 

-  In  TX  -  aXCL,  (B1 ) 


where  Tx  is  the  optical  transmittance  at  wavelength  A  (pm),  ax  is  the  mass 
extinction  coefficient  (m2  g~1),  C  is  the  aerosol  mass  concentration 
(g  »-3),  and  L  is  the  optical  path  length  (a).  It  is  straightforward  to 
show81  that 


aA  °  3Qx/20pP, 


(B2) 


where  the  Hie  theory83  gives  the  value  of  Qx  that  depends  on  the  particle 
diameter  Dp  (p»)  compared  to  wavelength  A,  and  p  is  the  particle  mass  den¬ 
sity  (  g  cm"3).  In  the  geometric  regime  (Dy»  A),  Qx-*  2.0  as  Op 
increases.  Thus. 


ax-*  3/Dpp.  (Dp»  A)  (B3) 


For  example.  Figure  81  shows  values  of  ax  vs  Dp  at  the  He:Ne  laser 
wavelength  (A  =  0.63  pm)  calculated  for  spherical  droplets  of  water  (dashed 
curve83)  and  several  phthalates  (dibutyl,  diethyl,  dimethyl,  dioctyl-  solid 
curves);  the  latter  are  commonly  used  in  aerosol  chamber  testing. BA 

The  figure  shows  that  when  Op  »  0.63  pm,  the  Hie-calculated 
extinction  coefficient  curves  merge  into  a  tail  whose  function  is  given  by 
Equation  (S3)  and  whose  thickness  depends  mainly  on  the  differences  in  mass 
densities  (p)  of  the  liquids  represented.  The  density  differences  are  not 
obvious  because  of  the  logarithmic  scale  used  for  the  ordinate  of  Figure  B1. 
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EXTINCTION  coefficient. 


Figure  81-  Hie-Calculated  Extinction  Coefficient  vs  Spherical  Droplet 
Qiaseter  for  the  He:Ke  Wavelength  (A  =  0.63  ths  Dashed 

Curve  Is  for  Water  Droplets,  While  the  Solid  Curves  Are  for 
Droplets  of  Phthalates  (Dibutyl,  Diethyl,  Dimethyl  and  Oioctyl). 
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A  great  many  liquid  and  solid  substances  exist  for  which,  when 
they  are  dispersed  as  spherical  particles  at  wavelengths  where  their 
absorption  is  negligible,  remarkably  good  first  approximations  of  a\  can  be 
made  using  Equation  (B3).  Hence  this  equation  is  useful  for  first  calcula¬ 
tions  involving  such  aerosols  when  their  optical  properties  are  not  known 
at  some  wavelength  k,  but  A  common  example  is  that  of  water  fog 

droplets  in  the  visible  wavelengths.  Often,  visible-wavelength  observa¬ 
tions  can  be  combined  with  IR  measurements  to  yield  additional  information 
about  an  aerosol,  e.g.,  mean  particle  diameter.85  Combining  Equations  (B1) 
and  (B5)  gives 


-  In  TX  =  3CL/D^p,  (DM»  X)  (B4) 

and  it  is  straightforward  to  show  that 

-  In  TX  =  Z[(Dp2  NL  x  10"6),  (D^»  k)  (B5) 


where  N  is  the  number  of  spherical  particles  per  cubic  centimeter. 

Equation  (B5)  can  be  extended  to  some  rather  practical  applications.  For 
example,  Middleton86  estimates  meteorological  range  as  that  for  which  a 
high-contrast  target  is  seen  with  2%  transmittance,  i.e.,  Tx  =  0.02  and 
-  In  Tx  =  3.912,  so  that  from  Equation  (B5j 

N  =  (2.49  X  106}/D^  L.  (Dm»A)  (B6) 

As  an  example,  a  typical  developing  water  fog  is  comprised  of 
droplets  having  mean  diameters  in  the  3-5  pm  range.  If  an  observer  could 
barely  see  a  high-contrast  target  at  L  =  200  m  under  such  conditions,  from 
Equation  B6  the  fog  should  contain  approximately  500-1400  droplets  cm’3. 
Perhaps  a  more  useful  application  would  be  to  calculate  mean  values  of 
from  Equation  (6)  using  observed  ranges  (L)  and  droplet  populations  (N) 
measured  by  independent  means,  as  a  water  fog  aged  and  dissipated. 
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